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Abstract. In Causal Perturbation Theory the process of renormalization is 
precisely equivalent to the extension of time ordered distributions to coincident 
points. This is achieved by a modified Taylor subtraction on the corresponding 
test functions. I show that the pullback of this operation to the distributions 
yields expressions known from Differential Renormalization. The subtraction 
is equivalent to BPHZ subtraction in momentum space. Some examples from 
Euclidean scalar field theory in flat and curved spacetime will be presented. 



Introduction 

Calculations in perturbative QFT are performed primarily in momentum space. 
The computation of a given contribution to the iS-Matrix is done by writing down 
Feynman rules and applying a certain choice of renormalization scheme to the 
resulting expression. For a reasonable renormalization scheme it should be proved 
to work to all orders and so produce a finite S-Matrix, as for example in the case 
of BPHZ renormalization. 

But today we consider the principle of locality to be of special importance, 
and hence a local formulation of perturbation theory should exist. Indeed this was 



elaborated upon by Epstein and Glaser |EG73| following earlier ideas of Bogoliubov 
|BS76 |. Their approach is called Causal Perturbation Theory (CPT). Based on a 
set of axioms they constructed the 5-Matrix as a formal power series inductively. 
The process of renormalization occurs only once in every step. All lower order 
contributions are already renormalized. This corresponds to the determination 
of all divergent subgraphs in the traditional approach and simplifies the proof of 
the construction to all orders. The main concept on which CPT is based is its 
formulation completely in configuration space. During the s eventie s there have 
not been many applications of it except in the works [ BS75 , DM75]. This may 
be due to the fact that Epstein and Glaser used rigorous functional analysis, so 
renormalization is defined by an appropriate subtraction on test functions, whereas 
physicists are used to working with distributions i n an in tegral kernel representation. 

Later, Scharf et al. applied CPT to QED (see |Sch95|| and references t herein) and 
to non Abelian gauge theories HDHKS94a| , [DHKS94b| , |E)HKS95| , |DHS95[| . But their 
perturbative calculations are still performed in momentum space. In my opinion 
CPT has much potential for further applications to computation of diagrams in 
configuration space. 

On the other hand, a renormalization scheme called Differential Renormalization 
|FJL92 , SZ93 | has gained attention. It is appreciated for its simplicity especially 
since no regularization procedure is needed. Differential Renormalization works in 
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configuration space on a large number of examples, but a proof for all orders is 
missing. 

In the present paper I will show that the integral kernel representation of the 
Epstein-Glaser renormalization procedure exactly yields Differential Renormaliza- 
tion. It will lead to a simple formula for the computation of diagrams. I will 
apply it to some examples from Euclidean scalar field theory and use the results 
for renormalization group computations. 

CPT mainly relies on the principles of causality, translation invariance and the 
singularity structure of the Feynman propagator. Brunetti and Fredenhagen [ BF97| 



implemented CPT on a globally hyperbolic spacetime by giving a local generaliza- 
tion of translation invariance. Here the local causality structure is preserved and 
the Feynman propagator is known to have Hadamard form. I will show how the 
corresponding distributions can be renormalized in the Euclidean case. This is 
achieved by an appropriate translation of their representations from flat spacetime 
to curved spacetime. 

1. The extension of distributions 



Following |Sto93, Fre96| it turns out that renormalization in CPT actually is 



an extension of distributions from the subspace of test functions whose support 
does not contain the origin to the space of all test functions. To treat the most 
general solution of that problem we will be concerned with the space of distributions 
V(R n ), the dual of V(R n ), the space of test functions with compact support. Let 
a = {ai, . . . , a n } € N™ be a multi-index, we set \a\ — Yh=i a i an d — n™=i a *'- 

Q\a\ 

d a = (1) 

d Xl ai ■ ■ ■ dx n a " {> 

is a partial differential operator of order \a\. 

Remark. Note that all operations on distributions like differentiation and trans- 
formations of their arguments are defined by the corresponding operations on test 
functions. 

This fact is referred to as "in the sense of distributions". Writing 

T((p) = f d n xT(x)<p(x), TeP'(M"), ^D(l"), (2) 

we call T(x) integral kernel of T. Let X>(R™ \ {0}) = {ip £ X>(R™)|0 g supp(y>)} 
denote the subspace of test functions whose support does not contain the origin 
and V'{W l \ {0}) its dual]] Now we state the 

Problem. Given a distribution °T £ 2?'(R n \ {0}), how can we construct an ex- 
tension T £ D'(R"), such that °T(tp) = T{tp) for ip £ D(R n \ {0}) ? 

The solution of this problem requires the introduction of a quantity that measures 



the singularity of the distribution at the origin |Ste71|. 



Definition 1. A distribution T £ V(R n ) has scaling degree s at x = 0, if 



inf{s' S R|A s 'T(Ax) ^ in the sense of distributions}. (3) 



1 The existence of the extension is guar anteed by the Hahn-Banach theorem. A solution for 
homogenous distributions can be found in ]Hor9c| ][Chap. 111.2] 
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Let s be denoted by scaldeg(T), and define singord(T) := [s] — n, the singular 
order 

The definition also holds if T G V'(R n \ {0}). Take the ^-distribution as an 

Example 1. For 5 G V'{R n ) one has: S(Xx) = \\\- n S(x). The scaling degree of 5 
is n, the singular order is zero. 

The scaling degree of some special compositions of distributions can be computed 
quite easily. We state 

Proposition 1. Let T G V'(R n ) or V'(R n \ {0}), scaldeg(T) = s and /3 be a 
multi-index. 
I. scaldeg(^T) =s-\/3\. 
II. scaldeg^T) = s + \(3\. 

III. scaldeg(w) < 0, scaldeg(wT) < s,w G V(R n ). 

IV. scaldeg(Ti T 2 ) =sx + s 2 , if scaldeg(Ti) = s h % = 1,2. 

The proof is skipped, we only note that all statements follow directly from the 
translation of the words "in the sense of distributions" and the use of the Banach- 
Steinhaus theorem (principle of uniform boundedness, applied to distributions) on 
point [II. 

Example 2. The scaling degree of 5^ G V'(R n ) is |a| + n. The singular order is 

14 

The solution of the problem depends on the sign of the singular order. Let us 
consider the simple case first. 

Theorem 2. Let a T G T>'(R n \ {0}) with scaling degree s < n. Then there exists a 
unique T G V'(R n ) with scaling degree s and T((p) = °T(ip) for all ip G X>(R n \{0}). 



The proof can be found in | BF97 |. If the scaling degree is not smaller than the 
space dimension, the singular order w is zero or positive. In that case theorem || 
guarantees a unique extension on test functions that vanish at the origin up to order 
ui. Thus a general extension can be defined after performing a projection into that 
subspace. This is achieved by a kind of modified Taylor subtraction, called the 

W- operation. Let T> u (R n ) be the subspace of test functions vanishing up to order 
uj at 0. Define 

W {ul , w) : V{R n ) i ► V"(R n ), if ^ W (u) . w) <p, 
(W (uiw) <p) (x) = f{x) - w(x) J ( d °B (0) ' ^ 

| a| <cj 

with w G V(R n ),w(0) ^ 0. 

The action of Wr u - W ^ on ip can be written as 

{W {U] ; W) <P) (X) = ^WW- 

\/3\=u+l 

with ipp G T>(R n ). It has the nice property 

W^:w)W(p = wW^-^ip. (6) 



[s] is the largest integer that is smaller than or equal to s. 
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With (d a x~<)(0) = j\SZ it follows for | 7 | < oj: 

W (u ., w) wx* = wW^xi = 0. (7) 
Now we can discuss the general case. 

Theorem 3. Let °T e D'(R n \ {0}) with scaling degree s > n. Given w G V(R n ) 
with w(0) 7^ 0, a multi-index a, \a\ < uj and constants C a G C, then there is exactly 
one distribution T' G T>'(R n ) with scaling degree s and following properties: 

I. (T',<p) = (°T,<p) V^e2?(K"\{0}), 
II. {T',wx a ) = C a . 

T' is given by: 

(TM = (T,W { ^) + E ^ ( da l) (0)- (8) 

|a|<w 

ifere T is the unique extension by theorem [|, W( u;w ) is given by (Qj and a; is the 
singular order of T. 



The proof can be found in JBF97|. We see that in the case of non-negative 
singular order, the extension is not unique. It is fixed by a finite set of complex 
numbers C a . Let us look at the next 

Example 3. The n-th power of the scalar Feynman propagator (iAp) n (x) = 
e(x°)A+ n (x) + 8(~x°)A + n {-x) is a distribution on R 4 \ {0}. We compute the 
scaling degree of A + ". 



A+ n (Aa;) = (2t 



A" 



n y3n f TT d3 Pi c T^-A-iu t , i Xx +iviXxi 

TT d P' P £r = i(-V(^") 2 +p» 2 ^ n +'P^) 

A i 2 A /(Am) 2 + Pi 2 



f \ 2^(Am) 2 +p 

= A- 2 "A+™(a;,Am) 

Here D + denotes the massless scalar two point function. Hence the scaling degree 
is 2n. The application of the VF-operation with u> — 2n — 4 yields the extension 
to all test functions. The computation can be done similarly for the Euclidean 
propagator. 

We turn to another example that seems to have caused some confusion in classical 
physics (see e.g. [FLS64]). 

Example 4 (The self energy of the electron). In electrostatics the electric poten- 
tial of an electron at the origin is given by the Green's function of the Laplace 
equation in 3 dimensions. 

Acj) = -4tt P = ineS => tf> = -- e VCR 3 ). 

r 

The electric field is 

E = -v</» = — ,e V'{R 3 ). 
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Since sing supp(E) = {0} it follows: 



E 2 = ^ G P'(M 3 \{0}). 



The singular order is 1. Hence there is an extension to all test functions by the W- 
operation. We can define the energy density U = E 2 as the following distribution: 

(U, f) := (E 2 , W(i ;w -)(f) . (9) 

An electron at rest has self energy E = 1/ (Air) (U,l). The choice of ip = 1 is possible 
due to sufficient convergence at long range. The same holds for w in (||). 

E = ^(U,W {1 . A) 1)+C = C° 1 

as all C a , \a\ — 1 vanish. We can determine C° by the requirement that the mass 
of the electron is purely electromagnetic, i.e. 

E = mc 2 . 

In the following I will suppress the distinction between °T and T in the case 
that the scaling degree is smaller than the space dimension. This should lead to no 
confusion since the extension is unique in that case. 

1.1. The integral kernel representation. To use standard vocabulary we will 
call the extended distribution in theorem |^ a renormalization. Next we will work 
out its integral kernel. If we set all C a to zero we have 

Definition 2. Let T G V"' (R n ) with singord(T) = uo. The integral kernel T R{uJ . w) 
of its extension is given by 

(T R(ul . w) ,<p}:=(T,W iu;w) <p). (10) 

Furthermore, we consider a family of distributions T t that depend continuously 
on a real parameter t. If K is a real compact interval then J dt (T t , if) exists as a 
Riemannian integral. We define 



(11) 



in the sense of distributions. Now we have 

Proposition 4. The integral kernel (Tm u . w \w\ is given by: 

(T MW )(,) = (-r + V + i) £ ^/"'dt^^r (12) 

Proof. The Taylor expansion of if at the origin is: 

*)=E^ )( ° ) + (tlJ + 1) £ i I d*(l-TOV)(ta). (13) 
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Hence {W[ UJ] i)fi){x) is the Taylor rest term of order ui+l. Writing ( |l2| ) as (Tmu^w) 
f dt (T R ^. w )w) t we find using @: 

Q dt(T R(w . w) w) t , i p S j = (u} + l) J di(l-t)"x 

x E ^r( a >( a: )(0'V)(te) 



= (T, W^. w )W(p) 

= (Tr( u - w )W, if) , 

where we used equation (p|). □ 



Note that the differential operator in (|12|) is a weak derivative. If w had no zeros 
(=> w ^ 2?(R' 1 )), the integral kernel would result from a simple division. Later we 
will encounter a well known example of this. 

To achieve a similar representation for the whole distribution we put a restriction 
on the test function w in (g). Let w(0) = 1 and (d a w)(0) = 0, for | a | < u>. 

Remark. This is no real loss of generality since for a given v € D(M n ), w(0) 7^ 0, 
the projection 

v ^ w = v ^j( da ~) (°) e V ( Rn ) ( 14 ) 

yields (<^w)(0) = for | 7 | < u. 

Then (1 — w) vanishes up to order w at 0: 

W {uJ . w) (l - w)if = (1 - w)ip, (15) 
((1 - w)T R(bJ . w) ,ip) = ((1 - w)T,ip) . (16) 

Now the integral kernel is given by 

Lemma 5. With the above restrictions on w <G T>(R. n ), the renormalized distribu- 
tion Tr( w;u ,) has the following integral kernel: 



T^)(aO = (-)> + !) £ dP i\-f At l^T T {l) W {l 

\P\=LO + l 



(i-fr /x\ fx 



+ / dt 



(1-*) 



(f)(-»)(f) 



(17) 



Proof. The first term of of ( |l7| ) is the integral kernel of wT R ^. w y A simple compu- 
tation shows that the second term smeared with <p yields (T, (1 — w)ip) . Equation 
161) completes the proof. □ 
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Y^d^xA \[ = ^d Ml x^ + ... + d flm x^)\ (18) 

101=1 / 

/ 'x» 

with x — 




e R 4 " 1 . 



Here \ . . . \ denotes the ordering of differential operators to the left of the coordi- 
nates. 

To perform some computations with formula ( |l7| ) it would be desirable to aban- 
don the requirement of w being a test function. Let w m G 2?(R") be a se- 
quence with lim m ^ooW m =: w G V(R n ). If lim m _,oo (Tk(w;w m ) j <£>) G C exists 
G P(R n ), we will allow w to be used in the renormalization procedure. Let 
us consider T 6 X> w (M n ) with singular order cj and singsupp(T) = {0}. Choose 
w(x) = 6»(1/A/ - |x|) =: 6<(x),M G R, M > 0, where | • | denotes the Euclidean 
norm. Since singsupp(T) n singsupp(0 < ) = and 9 < has com pac t support, the 
pointwise product 6<T G £ W '(R") C P W '(R") exists.| Applying @ yields 



i/3i~r + i p! Ji 1 yi/ as) 



For an arbitrary choice of w, a scale M has to be introduced for dimensional rea- 
sons. This allows writing w{Mx) as a function of a dimensionless argument. The 
dependence of the counterterms on the scale can easily be computed: 

M-^ (T R{uMMx)) ,</?)= £ Ba ( S{a) > V>) > (2°) 

\a\ <lo 

with B a = (T, Mx»(d u w)(Mx)x a ) . (21) 



For w = 8 < we get 



B- = U!l/r,5(i-|x|]y ). (22) 



1.2. Momentum space and BPHZ renormalization. Since the Fourier trans- 
formation is a map 2? ^ <S we have to restrict our distributional space to S' C 
We remind the reader that the Fourier transformation is defined in the sense of 
distributions, i.e. T(ip) := T(<p). 

First we start with the definition of the moments of a test function: 

K a (V0:= J d" x x a i/j(x), fle5(l"). (23) 
3 £ = C°° and £' is the space of distributions with compact support. 

4 Let S be the space of C°° functions of rapid decrease and S' its dual. We use the convention 



■$(p) = J d n xip(x) ( 



IpX 
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Let ^(R™) := {tp E S(M. n ), K a (i/j) = 0, \a\ < oj} be the subspace of test function 
with vanishing moments up to order ui, then it follows: tp E 6>" => ip E S u . Choosing 
w as in Lemma || we have: 

K a (w) = for < |a| < u, K°(w) = (2tt)", (24) 

and by a simple computation: 

K^(d a w) = (-) l7l 7 !^(27r)™, 7<a. (25) 

The Fourier transformation of Wip is: 

(w^, w)V y (p) = m - e i - z ^(^ a )(-p) (s {a \v) (26) 

= * {P) ~T^ E ^W(-p)- (27) 

1 ' |ot|<ui 

Using ( p5| ) we get: 

^ 7 ((^ ;ii ,)V) V ) =0. M<w, (28) 
so W( u . w y actually is a projector S i— ► 5 W . With 

(Tr( U ;w),¥>) = (Tr{u- w ),^ = (f, {W^. w) (fY^ (29) 
the integral kernel is given by 

tS^o w = f (*) - E 5 ( aQf ^) (°) ■ ( 3 °) 

|a|=0 

This subtraction should be understood in the sense of distributions, i.e. the sub- 
traction on the test functions has to occur before smearing out. It looks similar 
to BPHZ renormalization which is a Taylor subtraction at arbitrary momentum q. 
We will compute the corresponding w. Let 

n(k) ■= f(k) - J2 fc^( 9 «f )(q) (31) 

denote the BPHZ renormalized distribution in momentum space. Using the Taylor 
rest expression similar to ( fL3| ) and performing Fourier transformation we get: 

■)'-' ' 

|/3|=w+l 

and comparing to ( |l2"| ) , 

= T R ^. e iq^y (33) 

Here, the equivalence of the subtraction procedure in Epstein-Glaser and BPHZ 
renormalization can be seen explicitly. Moreover this cannot be achieved by the 
requirement w E S(R n ). 



Ti(.) = (-r+wi)e-- E d ^f Q d 4+^ T (!) e ^' (32) 
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2. Causal perturbation theory 



In the followi ng I wil l give a very brief summary of CPT. A complete description 
can be found in |]EG73| , |Sch95| . 



We start with the S-Matrix as a formal power series 

= 1 + 2_, — d 4 xi- ■ ■d 4 x n T n (x 1 , . . . ,x n )g(xx) ■ ■ ■ g(x n ). (34) 

n=l n ' 

It is an operator valued functional. The function g E 2?(R 4 ) plays the role of 
a coupling "constant". The T n are operator valued distributions in Fock space. 
They are called time ordered functions and involve free fields only. Epstein and 
Glaser stated six axioms from which the time ordered functions can be computed 
recursively. Here we cite only the most important one called the causality axiom: 

T n (xi, ...,x n ) = T k {xi, . . . , x k ) T n ^ k (x k+ i, . . . , x n ), (35) 

if all points x k +i , ■ • ■ , x n are not in the causal past of x±, . . . , xu . The inductive con- 
struction starts with T\ = C- ln t. We assume that T n i{xx, . . . , x n >) for all n! < n exist 
as a sum of products of a symmetric translation invariant numerical distribution 
and a Wick polynomial of fields. The scaling degree of the numerical distribution 
is known at coincident points. Now T n can be constructed up to the total diagonal 
xi = ■ ■ ■ = x n by the causality axiom. Wick's theorem ensures the required form. 
The numerical distributions have an extension to the diagonal which is the origin 
in R 4 ™ -4 because of translation invariance. The Wick polynomials are already de- 
fined as operator valued distributions on the whole space. Let me emphasize that 
the translation invariance of the numerical distributions plays a crucial role in the 
whole construction. 

3. Applications 

I will give some examples from 4 -theory in lowest order. The Lagrangian for 
the self interacting scalar field is 



Now T\ is given by the interaction term. 



£(x) = - .. dli cj>(x)d^(x):-— :4 2 (x):-- :^ 4 (x): . (36) 



T 1 (x)=~:4> 4 (x):. (37) 



Causality ( |35[ ) implies T 2 for non coincident points. 



T 2 (xi,x 2 ) 



:0 4 (x 1 )0 4 (x 2 ): + (38) 



' \Ti(x 2 )Ti(xi), ifx° 2 >xl 

A 2 
: (4!) 2 . 

+ 16A F {x 1 - x 2 ) :<l) 3 (x 1 )<j) 3 (x 2 ): + (39) 
+ 72(iA F ) 2 ( Xl -x 2 ) :(f) 2 {x 1 )(f) 2 {x 2 ):+ (40) 
+ 96 (»Af) 3 (si - x 2 ) \4>(xi)4>[x 2 ) : + (41) 
+ 24(iA F )\x 1 -x 2 )]. (42) 
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To give some explicit results for the numerical distributions we will now turn to 
their corresponding Euclidean counterparts. The singular support of these distri- 
butions is the origin only. Hence we can use (fill) for the renormalization. 



3.1. The massless theory. The Green's function of the Laplace equation in four 
dimensions is 

KM = (43) 

It has singular order —2. Consider the contribution to the two point resp. four 
point vertex function. 

3.1.1. The one-loop graph. The "fish" graph (Q) has singular order zero, so we get 

Df\ (x) = ^-fl^-L—i. (44) 
Since all extensions only differ by a <5-term we have 

M dM Df li {X) = 8^ S(X) 01 9 ^ = 27r2S{x) > (45) 

by (p2|). This also could have been seen by expressing t^jW = — 2^^) which is 
unique by theorem 0. 

3.1.2. The two-loop graph. The singular order of the "setting sun" (|il|) is 2. Hence 
we have 

D F 3 \*J (x) = -^—,\dn,d u d a x tl x u x' 7 T ^-r [ln(M 2 a; 2 ) + AfV - 4MVx* + 3 
lR (2tt) 6 4 M (x 2 ) 3 L 

(46) 

Terms two, three and four have singular order 0, 1 and 2. As they are zero outside 
the origin, they must be proportional to 8 and its first resp. second derivatives. As 
there is no Euclidean invariant combination of d and 5, term two has to be zero. 
Using ( p2| ) we get by comparison with M-^j on (|4^): 

d^d v d a x»x v x°-^ = ^AS(x), (47) 

d^d^x^x^-^-r = 4ir 2 6(x), (48) 

d fl d„d a x»x v x' r —^ s ee 0, (49) 



where A is the Laplacian. Now we will turn to 

3.2. The massive theory. The Green's function of the Euclidean Klein-Gordon 
equation is 

1 mKi ( mVx^j 

" (50) 

As Ki(x) <x 1/x for a; — * 0, the singular order of = —2. We will compute 
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3.2.1. The one-loop graph. With singord(Ai? ) = and (|19[) we get: 



1 Q J" 1 ' 

rd L ,X>*< 



32tt 4 M 



ft 2 ( mVx 2 ) — Kq mV? if; 



rn 2 1 

M 2 (a; 2 ) 5 



ft 2 



2 mv a; 



ft r 



(51) 



By writing x^ ■ ■ ■ = . . . which is unique we can express it as: 

ft (mVx 2 ) ft mV? 



2(2tt) 4 



Ai ft 2 (mVa?) - ft 2 (mVa?) 



1 



2M 2 x 2 



ft 2 



MS) 



ft - 



VM 



(52) 



This can be compared to the corresponding expression in | HL92 |. 



3.2.2. The two-loop graph. Applying formula ( |l9| ) leads to an integral J ds (s — 
const) 2 K \ (s) that is difficult to solve. But if we use the expansion 

m 3 K 1 3 (mr) 1 3m 2 / ( m 2 r 2 \ , , \ , „, ,„ , 
= ^ + 4^ (—J + ln ^ ) - V + fl ^ r )' ( 53 ) 

we can renormalize the first and second summand with singular order 2,0 respec- 
tively. With 



In 



(x 2 ) 



2\2 



M 



1 ln(i^)ln(M 2 ^ 2 ) 
4 " (x 2 ) 2 



(54) 



and the previous results from subsection |3jj we get: 

(A F 3 )C R (x) = * S^d^d^W-^ [HM 2 x 2 ) + M 2 x 2 + 3] + 



(27T)' 
1 

16 "^(^2)2 



+ —9,,,/ 



In 



4 2 

m x 
V 16M 2 



21n( 7 2 



ln(AfV) + i?(x" 2 ) L (55) 



where the subscript MR denotes the "minimal renormalization", i.e. every summand 
is renormalized with its singular order. 



3.3. The renormalization group. The parameter M plays the role of a renor- 
malization scale. It enters the theory by purely dimensional reasons as an argument 
of the "function" w. We will give the lowest order contributions to the 0, 7 and 
7,„ functions in the renormalization group. They can be read off by solving the 
renormalization group equation 



r„(:E] 



j . . . j ^ n j 







(56) 
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to order g 2 , where g = A/(167r 2 ). Using the expansions 



p(g,m,M) = M-gjL=J2Pn9 n , 

n=2 

dlnm 



r ym(g,m, M) = M 
l(g,m,M) 



dM 
M dlnZ 



n=2 



2 dM 



E 

n=2 



7n9 



(57) 
(58) 
(59) 



we find /?2 = 3 and 72 = f° r the massless theory. Here we had to add a term to 
the setting sun proportional to ([Is] ) to achieve 7 m = 0. 
In the massive theory we find 

™2 



02 = — 



M 2 



M 



(60) 



hence 02 — » 3 if M — > oo. This result also holds for the corresponding minimal 
renormalization. In that scheme we get 



72 =12' 



7m,2 



1 

2M 2 
3m 2 



In 



/ m 2 \ 
I 4M 2 ) 



Hi) 



5 

12' 



(61) 
(62) 



by using (|5 



3.4. Curved spacetime. Let M. be a globally hyperbolic manifold with a metric 
g. Wick polynomials were defined in ||BFK96| using techniques from microlocal 
analysis. Then CPT was implemented in | BF97 | for scalar 4 -theory. The Feynman 
propagator is known to have Hadamard structure | KW91 |: 



A* 

Af oc — h uln(2cr) + w, 

2a 



(63) 



where a, A, v, and w are smooth functions on the manifold and an appropriate 
ie regularization has to be chosen. By using a chart it can be seen to have the 
same scaling degree as in flat spacetime. The function a is half the square of the 
geodesic distance which is unique in every sufficiently small neighbourhood. Let 
g = det(g a b). The Van-Vleck-Morette determinant 

A ^' p ') = / 7 w „ det(-<7 n y (p,p')) (64) 

V9{p) g(p) 

fulfills the following differential equation: 

V a {Aa a ) =4A. (65) 
The vector index on a denotes the covariant derivative as usual. If we expand 



E 



w 



E 



w n a 



(66) 



in powers of cr, the coeffic ients (except wq) are determined by the Hadamard recur- 
sion relations, see | DB6C |. 
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Following | BF97| the causal construction of the second order 5-Matrix is the same 
as in flat spacetime. The only step left is to perform the extension to coincident 
points. 

Therefore we turn to the corresponding Euclidean distributions, so we ca n trans- 
fer our results from the previous examples.^] They may be compared with | Liis82 
Bun81|. Again we only use the minimal renormalization scheme. 



3.4.1. The one-loop graph. This is given by 



■ A/ 
\MR 



16tt 4 



(2af 



A I 



2 rr-i — '- + 2 — h 2vw \n(2a) + v 2 In (2a) + w 2 



2(7 



2d 



Only the first term needs renormalization and we get 



A 



(2af 



At 



1 Aln(2M 2 q) 
2 Vaa (2af ■ 



(67) 



(68) 



3.4.2. The two-loop graph. Here we use the expansion ( |66| ) to determine the terms 
that have to be renormalized. Then we have 



A P 3 



M 
MR 



1 



At 



A I 



(2tt) 6 \ (2a) 3 
At 



Av ln(2cr) 



(2a) 



A I 



Aw 



(2a) 



A I 



A» ln(2er) 
' 2a 



_ A2 V 2 ln(2a) A2 UW ln(2cr) A^w 2 

3— + 3 — - — '- + 6 — ^ — '- + 3— — 

2a 2a 2a 2a 



+ v 3 \n 3 (2a) + 3v 2 w \u z (2a) + 2vw z \n(a) + w :i 



(69) 



with 



1 °° 



w 



n=0 

The first term is found to be 

M 

1 



n=0 



A3 



(2af 



= -(V a V fc V c A2 - 3V a V h Ai c + 3V a Ai &c - A^ , abc )a a a b a' 



(70) 



ln(2M 2 cr) 



n 



(2af 



Similarly the second and third terms can be computed: 



Av \n(2a) 



(2aY 
Aw 



(2af 



I. , Aln(2M 2 a)ln 
= -(W a vo-v , a )a 

l, s a Aln(2M 2 (j) 

= -(V a w -w 0;a )a (2g)2 ■ 



' 2a ' 

. M 2 , 



(71) 



(72) 



(73) 



Without further knowledge about the Feynman propagator, the /3-function can be 
evaluated to lowest order, leading to the result /?2 = 3 in this renormalization 
scheme. The calculation requires the use of the identity V a = 27r 2 <5(p,p'). 



For the a calculus see e.g. [Ful8E] 
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Conclusions 

The elegant method of CPT is suited not only for the investigation of renormaliz- 
ability but also for the performance of perturbative computations. The subtraction 
procedure on the test functions can be pulled back to the distributions yielding Dif- 
ferential Renormalization. Therefore it is possible to work in the standard integral 
kernel representation. 

As the whole procedure is formulated in configuration space, it can be transfered 
to distributions on a manifold. This enables one to give compact expressions for the 
renormalization of quantum fields in curved spacetime, at least in the Euclidean 
case. 
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